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I. INTRODUCTION 



Particle dynamics in a strong static magnetic field implies a separation of scales that 
allows for a reduction, the so-called guiding-center reduction, which removes the fast time- 
scale from the dynamics 1 . The fast time-scale relies on the gyro-angle, which is the angle 
measuring the Larmor gyration. The principle is to remove the presence of the gyro-angle 
in the dynamics of the other coordinates, by building suitable coordinates order by order in 
a small parameter expansion, linked to the Larmor radius. 

Guiding-center theory is widely used in plasma physics, because many physical phenom- 
ena rely on time scales much slower than the gyro-frequency, and spatial scales much bigger 
than the Larmor radius. For instance, Gyrokinetics is a kinetic model of plasma dynamics 
induced by guiding-center theory, and it is a key model for the study of magnetically confined 
hot plasmas, where it has proved efficient for the study of plasma micro-turbulence 2 . 

The coordinate transformation is usually identified by Lie-transforming the phase-space 
Lagrangian. Here, we consider performing the reduction by Lie-transforming directly the 
equations of motion. Although it does not guarantee a Hamiltonian structure for the slow 
reduced coordinates, it brings out a simplified derivation, while dealing with the gyro-angle 
in a more intrinsic way. 

Indeed, the guiding-center reduction at first order was already known a long time ago, 
but an efficient method to go to arbitrary order in the Larmor radius expansion was more 
difficult to get. 

Alfven already studied the lowest order slow particle drifts across the magnetic field 
lines, when the fast gyro-motion is removed 3 . Later on, particle dynamics at first orders 
in the Larmor radius expansion was computed 4-6 , especially because it provided a more 
precise expression for the dynamics of a virtual particle, the so-called guiding-center, which 
has only a slow time-scale motion, and for an adiabatic invariant, the magnetic moment. 
Those results were obtained working on the equations of motion, expanding them order by 
order in the small parameter and removing the fast time-scale through an averaging process. 
Computations were very heavy : to reach part of the components involved in the second 
order reduction, the amount of algebraic manipulations was so huge that mistakes could 
not be avoided 6 , and one could not consider going further in the derivation without a more 
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systematic procedure. This was all the more necessary as a wide arbitrariness was beginning 
to appear, and it was not clear how to deal with it in a natural way as had been done at 
first order. This required to make the derivation rely on a clear mechanism explaining what 
this arbitrariness meant, and what it could be used for 6 ' 7 . 

As particle dynamics in a magnetic field is Hamiltonian, Hamiltonian systems were good 
candidates for the required systematic procedure, with their quite developed perturbation 
theory. The guiding-center reduction was tackled with canonical coordinates 8,9 , which are 
more common and have richer reduction tools. Such a description did not deal with the 
magnetic field, the physical momentum and general spatial coordinates, but with the vector 
potential, the canonical momentum and magnetic field-lines coordinates. It did not easily 
fit in with non-Hamiltonian approaches. Also the separation between the unperturbed and 
the perturbed dynamics was not satisfactory. 

A breakthrough came with Littlejohn's works 7,10 . He used not only a Hamiltonian descrip- 
tion, but a non canonical description, because the physical particle coordinates (p, q) are 
non-canonical. In addition to classical tools, such as the Darboux theorem and expansions of 
differential equations, he used (canonical) Lie-transforms 11,12 . This Hamiltonian procedure 
was systematic, more suited to the time-scale separation, and closer to non-Hamiltonian 
approaches. The efficiency of the method allowed to include as easily the case where an 
electric field is present. Still, the procedure remained heavy, especially because it included 
several changes of coordinates since it proceeded in two steps, working first on the Poisson 
bracket, and then on the Hamiltonian. 

Finally, Littlejohn condensed his method into a single non-canonical change of coordinate 
by Lie-transforming directly the phase-space Lagrangian, which contains both the Poisson 
bracket and the Hamiltonian 13 . It was an efficient procedure to go to higher order. Since 
then, guiding-center reductions almost often use Littlejohn's approach 1 , and such an effi- 
cient procedure is often considered to be a result of the Hamiltonian reduction tools, which 
explains why Northrop's approach working on the equations of motion was so difficult, if 
not impossible at order higher than two. 

Actually, for the removal of the fast time-scale from the dynamics, the key ingredient 
is just the Lie-transform. It is what makes the reduction process systematic, as already 
appears through Littlejohn's work 12 . A similar idea can be found in the work 14 , which 
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worked on the Vlasov operator and expanded it to identify the guiding-center transformation. 
The Hamiltonian character is useful, for instance to guarantee the reduced dynamics to be 
Hamiltonian, with all the associated conservation properties, but it is not essential. The 
reduction can be efficiently performed working directly on the equations of motion. 

The resulting mechanism is more straightforward than when working on the phase-space 
Lagrangian : the whole second order reduction is obtained in a direct computation, despite 
the complicated expressions of the second order transformation, and the algorithm is so 
simple that it can be very easily implemented in a computer to reach any order. 

It also comprises an interesting approach of the guiding-center arbitrariness. Indeed, the 
guiding-center reduction is not unique, as is again emphasized by some recent results 15 ' 16 . 
Such an arbitrariness caused troubles in initial works dealing with the equations of motion, 
where it was unclear what the natural choices were to be at second order 6 . Hamiltonian 
approaches did not really address the arbitrariness in the averaging procedure, because 
they average not just the dynamics of the four reduced coordinates, but the whole reduced 
Lagrangian, with its seven components, including the Hamiltonian and the Poisson bracket. 
Now, when Lie-transforming the equation of motion, the arbitrariness is distinct and the 
natural choice consists in setting to zero all arbitrary terms in the transformation generator, 
exactly as it was done at first order to define the guiding-center position. 

The resulting guiding-center reduction is unique and minimal. The requirements are just 
tailored to reduce the dynamical dimension, which means just to extract a slow dynamics 
for the four reduced coordinates, the guiding-center position and the parallel velocity or 
the pitch-angle (the angle between the magnetic field and the particle velocity). Then, 
the method is short from the requirements to the related equations to be solved. It just 
corresponds to making a near-identity change of coordinates, i.e. a Lie-transform of the 
velocity vector field. Computations are straightforward and fully algorithmic, since the 
equations are solved just by an explicit inversion of the Larmor gyration generator. Last, 
the result is the minimal transformation that allows for the desired reduction. Especially, 
the vector field generating the transformation is a pure gyro-fluctuation. 

In the derivation, we avoid introducing the usual gyro-angle, with its issues about gyro- 
gauge 17 ~ 21 . Only the four reduced coordinates have to be changed, since only they have 
to be gyro-averaged. For the remaining coordinates one can deal with the initial physical 
coordinates. 
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The paper is organized as follows. In sect. 2, the goals of the reduction are introduced 
and written as equations for the reduced coordinates. It is shown that these equations 
can be solved explicitly at any order in the gyro-angle, and the iteration mechanism of 
the reduction is identified. In sect. 3, it is shown how the reduction is computed and 
the resulting algorithm is given. In sect. 4, the result is written to full second order in 
the Larmor radius, obtained in a straightforward computation, a comparison is done with 
previous guiding-center reductions, and some insights are brought into the general structure 
of guiding-center results. 

In sect. 5, a few possible extensions of the method are studied, such as the presence 
of an electric field or non minimal guiding-center reductions. For instance, averaging also 
the dynamics of 9 is considered, as well as including the magnetic moment in the reduced 
coordinates or providing the reduced dynamics with a Hamiltonian structure. This will 
clarify the respective advantages of Lie-transforming the equations of motion and the phase- 
space Lagrangian. 



II. LIE-TRANSFORMING THE EQUATIONS OF MOTION 

The dynamical system is simply a charged particle with position q, momentum p, mass 
m and charge e, under the influence of a static inhomogeneous magnetic field B. The motion 
is given by the Lorentz force 

1 TO ' 

p = p x eB . 

For the sake of clarity, we consider no electric field. The presence of an electric field 
satisfying the guiding-center ordering would not change the method (see sect. 5). 

When the magnetic field is strong, the motion implies a separation of time-scales. This 
is best seen by choosing convenient coordinates for the momentum space 

P ■= IIpII , 

(f := arccos ( g ) , (1) 



c • iip_ 



Px 
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where b := is the unit vector of the magnetic field, and p^ := p — (p • b)b is the 
orthogonal projection of the momentum onto the plane perpendicular to the magnetic field. 
The coordinate p is the norm of the momentum, ip is the so-called pitch-angle, i.e. the 
angle between the velocity and the magnetic field. Following Littlejohn's notations 7 ' 10,13 , 
the vector c is the unit vector of the perpendicular velocity. 

Differentiating formulae (1), the equations of motion are found as 

q = 

p = 0, 

<p = -P . Vb-c, 

c = -—a - £- • Vb • (cb + aa cot u>) , 

mm V 1 TJ i 

where p is now a shorthand for p(bcos</? + csiny?), B is the norm of the magnetic field 
and a := b x c is the unit vector of the Larmor radius, following Littlejohn's notations, so 
that (a, b, c) is a (rotating) right-handed orthonormal frame. The computation of c is rather 
involved, but it can be avoided by using formula (12). 

In the case of a strong magnetic field, the only fast term, the Larmor frequency uil :— — 
concerns only one coordinate, c, the direction of the vector p^ in the 2-dimensional plane 
perpendicular to the magnetic field. This corresponds to an angle, the so-called gyro-angle, 
and measures the Larmor gyration of the particle momentum around the magnetic field. 

To get a true scalar angle instead of the vector c, one chooses at each point q in space a 
direction which will be considered as the reference axis e^q) e B ± (q). Then, the gyro-angle 
9 is defined from the oriented angle between the chosen reference axis ei(q) and the vector 
c : 

c = — sin 9ei — cos 9e 2 , (2) 
and the equation of motion for 9 is 



9 = — + cot if — ■ Vb ■ a + — ■ Vei ■ e 2 



with e 2 := b x ei is the unit vector such that (b, ei, e 2 ) is a fixed right-handed orthonormal 
frame. 
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To emphasize the fast term of the dynamics in strong B, the equations of motion can be 
expanded in B~ x . Writing z := (q,p, <p, 9), we get 



eB 



o 

1 I P sm <£ 

J I + eB 



b cot ip+c 


— (b cot <£+<:)• Vb-c 
(bcot <£+c)-(cot </pVb-a+Vei-e2) 



= Z_i + Z , 



with z_i := 



z := 




p. 

m 


-£-Vb-c 
m 



\ 



(3) 



Y ^■(cotipVb-a+Vei-e 2 ) / 



where the indices correspond to the order in B 1 . The term of order B 1 has dimension of 
— , whereas all the terms of order B° have dimension of — • V, if we use the correspondence 
between the vector field q and the differential operator q • V. 

The ordering parameter, that is the ratio of two following terms in the perturbation 
expansion, is then e = = ps e '°^ V, which is the dimensionless quantity r £ V, with r L := 



ps "g ip the Larmor radius. This parameter is well-known as the magnetic inhomogeneity at 
the scale of the Larmor radius, because in r^V, the gradient always acts on the magnetic 
field, which is the only local property of the configuration space. 

The ordering parameter is sometimes considered just as the Larmor radius r L , or as 
B -1 , in agreement with an expansion in strong magnetic field. Also, it is often considered 
symbolically as e _1 , because it is an equivalent expansion but it avoids to deal with a 
space-dependent parameter involving B ; this is symbolic, because e -1 is not dimensionless. 
Notice that the ordering parameter r £ V is not just a scalar but an operator ; it has only a 
kind of dimensional meaning : the terms of order e 2 may not be proportional to (r^V) 2 , but 
they will have the dimension of ^V 2 . This point will be illustrated by the results of sect. 4. 



The goal is to isolate the dynamics of the slow variables from the fast variable, i.e. 
to perform a change of coordinates r : z — > z such that the dynamics of the remaining 
coordinates (q, (p,p) does not depend on 9. This is already obtained for p, so one only has to 
change coordinates on q and if. The coordinate transformation r transfers to functions by 
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duality through the "push-forward" operator T 1 , defined by the scalar invariance property 1 : 

(T- 1 /)(z) = /(r- 1 z). 

At the lowest order in the Larmor radius e~ x , the requirements are trivially satisfied. 
So, the transformation can be near-identity. It can be written as the exponential of a Lie- 
transform z = e~ x z, with — X a vector field, generator of the diffeomorphism r, which 
satisfies 

X p> 6> = , 

since only the coordinates q and ip need to be changed. The index notation is used to indicate 
the components, e.g. denotes the component ip of the vector X, and X P) ^ denotes a 2 
dimensional vector, whose coordinates are the components p and 9 of X. Through the 
transformation, the equations of motion become 

z — > z = e L z , 

with L := L x the Lie-transform along the vector field X. 

Now, the goal is that the equations of motion for q and f do not depend on 9, which 
means that all their non-zero Fourier components (i.e. purely oscillatory terms) are zero : 

osc(z) q ^ = 0, 

where following Littlejohn's notations, osc = 1 — avg is the projector onto gyro-fluctuations, 
with avg the complementary projector onto gyro-averages : 

avg(/) = I PdB f 
Jo 

for any function /. 

Last, L is expanded in series in the small parameter e 

= osc (z) = osc (e Ll+L2+ -z) , (4) 

where L is of order e\ 
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Equation (4) is the equation to be solved for the change of coordinates L n+i . Expanding 
it in series in e, we get an equation for each order : 



= osc(z_i) qj¥ , , 

= osc (Uz_i + z ) qjV , 

= osc ( L 2 z_i + ^-z-i + Liz ) , 



(5) 



At each order n ^ 1 in e, the highest order unknown L n+ i is involved only in one term. 
Isolating it, the equation writes 

-OSc(L n+1 (z_ 1 ) q ^) = (N n ) q)V , (6) 

where N n is a shorthand for all the terms of equation (4) that are of order n and that do 
not include the unknown L n+1 , e.g. N x := osc (^-^-l + Liz j 
The operator to be inverted is 



-OSc(L n+ i(z_i)) q)¥ , 

= -OSC ((X n+ i)j<9i(z_i) qiV - (z_i)j5i(X n+ i) q)¥; ) 

= OSC ((z_i) e <%(X n+ i) q ^) 

= ^^(X n+ i)q )¥ , , (7) 

where, in the first equality we used the usual formula (18) for the Lie-transform of a vector 
field, in the second equality we used (z_i) qj¥ , = 0, and in the third equality we used that 
osc do is just do, because do takes its values in the gyro-fluctuations, i.e. the non-zero Fourier 
component in the variable 9. 

The operator u^do is the generator of Larmor gyration, with the Larmor frequency as a 
coefficient. Equation (6) has a solution because the right-hand side (N„) q ^ is in the range 
of the operator do, since it is a pure gyro-fluctuation. The operator is easily inverted as 

(X„ + i) q>v = avg ((X n+1 ) q , J +±jdB (N n ) q ,„ , (8) 

where avg(X n+ i) q( ^ is a free element in the kernel of do, that is a free gyro-averaged 
function. We defined f d6 (N n ) q itp as the primitive of (N n ) qiV with zero gyro-average, i.e. 



J dO (N n ) q)¥ , := osc (AT) with Af any primitive of (N n ) qj , 



As a result, equation (5) can be solved to arbitrary order in the small parameter e through 
formula (8). At each order, one only has to write formula (5) to the n-th order, group into 
(N n ) qj¥ , all the terms that depend only on quantities that are already explicitly known, 
expand the result, and last invert the Larmor gyration operator through formula (8). 

This is close to Kruskal's averaging procedure 4 , whose principle was to expand the func- 
tion defining the motion and average this expansion. Here, we focus on the vector field 
defining the motion and Lie-transform it so as to make it independent of the fast coordi- 
nate. Such a procedure actually does not rest on averaging methods so much ; it could rather 
be related to normal form methods 22 , which aim at giving a simplified form to a vector field, 
when studying the local dynamics, the stability and possible bifurcations of an equilibrium 
point. 

In the guiding-center transformation, at each order, the fluctuating part of X n+i given by 
formula (8) is necessary and sufficient to solve equation (4). The averaged part is completely 
free, but also completely useless to solve equation (4). So, a natural choice is to put it to 
zero 

avg (X n+1 ) = . (9) 

This makes the transformation unique : the result is the minimal guiding-center reduc- 
tion, i.e. the only transformation which gives the desired requirements and whose only 
non-zero components are the fluctuating part of (X n+1 ) q (/ ,. 

In this guiding-center reduction, the transformation affects only the coordinates (q, <p). 
The transformation generator has no gyro-angle component. Now, the corresponding co- 
ordinate 9 is not intrinsic. It implies a non-trivial gauge fixing, whose global existence 
can fail 18,21 . It is not purely physical : in all physical results, what appears is the physical 
quantity c, as is clear in the literature 1 ' 7 ' 13 . The use of 9 is a detour and should be avoidable. 

So, we remain with the variable c with its initial definition : the unit vector of the 
perpendicular velocity 

c : = -Ex = Ed£*h . (10) 
With this variable, the coordinate space is constrained : the gyro-angle c is not inde- 
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pendent of the spatial position, since c e b , and b depends on q. When the q coordinate 
is changed, the c coordinate cannot be kept unchanged, otherwise it may get out of b -1 . 
Differentiating relation (10) with respect to q, we find 

Vc = -Vb • (cb + aacoty?) . (11) 

This formula can be obtained more easily by noticing that in the change of coordinates 
(q, P) — > (<l,P, c ), we have 

— sin pVp = V cos p = Vb • E = Vb • c sin p 
^ Vy? = -Vb • c 

=>• = V (j^j = Vc sin ip + Vb cos p + Vy?(— b simp + c cos p) 

= (Vc + Vb • cb) sin p + Vb ■ (1 — cc) cos p 
Vc = -Vb • (cb + aacoty?) . (12) 



in which V means differentiation with respect to q while keeping p constant, and we used 

= 0. 



that 1 = aa + bb + cc and that b is a unit vector, which implies Vb • b = V y-j 

The action of V on the vector c must be taken into account through formula (11) when 
computing the action of z q or (X n+1 ) q , but also when defining the components of z : writing 
/ = z • d z f in coordinates z := (q,p, p, c) with the property (11) implies 



z = 



eB 





] _|_ P S1 " r 



eB 



= Z_i + Z , 



with z_! :- 



z :-- 




b cot ip+c 


-(b cot <p+c)-Vb-c 




(13) 



Notice that the c component of z is zero, because the term — ^ • Vb • (cb + aa cot p) in the 
dynamics of c comes from z q • V/ with formula (11). 
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III. THE REDUCTION ALGORITHM 



In this section, it is shown how the computation proceeds for the minimal guiding-center 
transformation. To lowest order e^ 1 , the equation to solve (5) writes 

= osc(z_i) q>v . 

From the definition (3) of z_i, this is trivially verified, and it is actually a condition for 
the near-identity Lie-transform to isolate the fast time-scale, which is possible only because 
at lowest order, the motion concerns only the fast variable c. 

To order 0, the equation to solve (5) is 

= osc (L 1 z_ 1 + z ) qiV , • 

The solution for the first-order change of coordinates X x is given by equation (8) with 
the choice (9) 

( X iC = ^ J Moscow (14) 
The spatial component gives the lowest-order Larmor radius, which is often identified 
with the Larmor radius itself 



(Xx) q = ^Jd6 osc(z ) q 

= 7b j " osc (b cosy? + csin</?) 

= £ JdBcsm<p=e%£a = r L a, (15) 

where the action of osc and f d9 can be computed trivially by using formula (2). One can 
avoid this trick and make the computation with purely intrinsic operators using the results 
of 23 . 

On another hand, the component (p of equation (14) gives the expression of the first-order 
change for the coordinate ip 

(X 1 ) ip = ^Jd6 osc(z ) v 

= — ^ J d6 osc((b cos ip + csin <p) ■ Vb ■ c) 

= —jL J d6 (cos cp b ■ Vb • c + 2 sin y?Vb : a 2 ) 

= - E S £ (cot ¥ ?b-Vb-a-iVb:a 1 ) , (16) 
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with ax := — ^jr 2 and a 2 := S£z ^ ± the standard dyadic tensors of guiding-center works 1,7 . 
Formulae (15)-(16) agree with the usual fluctuating first order generator of the guiding- 
center reduction 1 ' 7,13 . 



To order one, formula (5) writes 



|_2 

= OSC L 2 Z_i + -^Z_i + Uz 



zo) 

As usual, it is already solved by formula (8) with condition (9) 

(X^i/d* obc (fz^ + Uzo)^. 



(17) 



All we have to do is to make the left-hand side of equation (17) explicit. This is completely 
algorithmic. The first step is to compute the Lie derivatives using the standard formula 



L n w = L Xn (w fc d fc ) 

= ({X^jdjWk - WjdjiXn^dk + (X n ),,w m [dj,d m ] 



(18) 



for any vector field w. Einstein convention is used and repeated indices are implicitly 
summed. The derivative operator <9 fc corresponding to the gyro-angle variable c is p- b x d p = 
— a • d C: since it is the generator of Larmor gyration 23 , and is equal to dg. 

In formula (18), the commutator of derivatives [dj,d m ] appears because the coordinate c 
is constrained, it is space-dependent, and the corresponding connection involves the pitch- 
angle. So, the following commutators are non-zero : 



[V, -a • d c ] 
[Vi,V,-] 



cot if (Vb • c) (a ■ d c ) , 

-(1 + cotV)(Vb- a) (a-a c ), 

(1 + cotV) f(Vib-c) (Vj-b-a) 



(Vj-b-c) (V<b-a) (a-d c ) 



Formula (18) is systematically applied to LiZo, then to l_iz_i, and last to (LiZ_i), 
which appear in the right-hand side of (17). Although uncomplicated, computations must 
be very orderly to remain tractable since the products and the Leibniz rule generate many 
terms from a single expression such as ^rz_i . 
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The second step is to perform the action of the gyro-integral J d6 osc. It acts only on 
the variables c and a, involved in expressions such as 

a • VB c • Vb • a = (aca) ijfc (Vfl Vb) ijk , 

upon which f d9 osc operates on the first tensor in the left-hand side by mixing the c and a. 

A way to perform the action of J dd osc is to compute its action on the basic tensors c, 
cc, ccc, etc. and to deduce its action on the other tensors by using cross products with the 
magnetic field (bx)^ = e^bfe ; this last operator is usually denoted by b^. For instance, 
(aca) can be written 

(aca)jjfc = bj m bkn (ccc) m j n . (19) 

Last, the action of J dO osc on the elementary tensors c® N can be computed by using the 
intrinsic calculus introduced in 23 , or by introducing a local fixed basis (ei,e 2 ) through the 
change of coordinate : 

c = — sin 9ei — cos 9e 2 , 
a := cos#e! — sin6>e 2 , 

as is standard in guiding-center reductions. With the fixed basis, the action of f d6 osc 
is trivial to compute : osc just cancels the zeroth Fourier component (gyro- average), and 
f dO is an easy integral. If iV is large, the computation by hand may be tedious but it 
remains trivial with a computer. Then one can come back to the initial basis (c, a) to get 
Jd6oac(d* N ). 

For instance, the lowest orders formulae are 

dOosc (c) = a , 

C3 -4- 3C 

d6osc (cc) = j , (20) 



/ 



J dOosc (ccc) = | [acc + cac + cca + 2aaa] , 
dOosc (cccc) = 5(accc + cacc + ccac + ccca) 

+ 3(caaa + acaa + aaca + aaac) 



/ 



These are all we need to get the second order reduction. Higher harmonics could be as easily 
dealt with. In practical computations, the relations given above can often be simplified by 
the symmetries of the tensor which c®^ is contracted with. 

14 



Notice that f dOosc is a linear operator which preserves the order in the fast variable c 
or a. It is a pseudo-inverse for dg. If we restrict the operator dg to the set of tensors that 
are harmonics of order n in the fast variable, it becomes just a matrix with finite dimension 
2™. If n = 2N + 1 is odd, the matrix is invertible, and is easily obtained (on a computer 
for instance) and gives the action of J dOosc. As regards even harmonics n = 2N, the 
kernel of dg is not zero, but it is a complementary space to the range of dg, as is obvious in 
Fourier series, for instance. So, after identifying the kernel and the range of this operator, 
one gets an invertible matrix by restricting the operator to its range. Then, the action of 
the operator fdOosc is just this invertible matrix on Range (cfo) and zero on ker(<%). It is 
actually a very efficient way to perform the action of the operator J dOosc. 

After using formulae (19) and (20) on each term of -yZ-i + Uzq, recombining the various 
terms and simplifying the result with formula 

b- VB = -BV-b, (21) 

coming from V-B = 0, one obtains the second order generator of the minimal guiding-center 
transformation X2, which is given by formulae (22)- (23) in the next section. 

The result is somehow intricate, especially for (X 2 ) ¥ >, but this was expected for a formula 
at second order, which even was not reached in usual derivations of guiding-center reductions. 
Actually, it emphasizes how efficient Lie-transforming the equation of motion is, since it 
provides by a straightforward computation such a complicated result. 

The explicit computation at order 2 illustrates the method to perform the derivation to 
arbitrary order. To each order in the Larmor radius expansion, the algorithm consists in 
writing down formula (8) and making all the terms explicit. Only two special operations 
are involved : the Lie-derivatives in the computation of (N n ) qi¥ ,, given by formula (18), and 
the gyro-integral f d6 osc computed by (20) ; these are automatic operations. 

If computations soon become tedious, it is only because the numbers of terms rapidly 
increases, as a result of products in Lie derivatives, of the Leibniz rule for the action of 
gradients, and of the gyro-integral operator. This is what generates the complicated formula 
(23) for (X 2 ) ¥ ,. But this is no trouble, since computations are completely algorithmic and 
very direct : formula (23) for (X 2 ) v , was obtained with no trouble by hand, and computations 
at higher order can be very easily done by computer-assisted symbolic calculus. 
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The fast growth of the number of terms with the order of expansion raises the question 
of the convergence of the series, but in perturbation expansions, convergence is not the first 
issue 12 and a factorial growth is quite standard. Even if the convergence is not guaranteed or 
when non- convergence is proven, the asymptotic expansion, with its truncated perturbation 
expansions, allows for a strong reduction of the effect of the fast time-scale in the dynamics 
of (q, </?), as emphasized in 4 . 

From this point of view, a simplified derivation, such as the minimal guiding-center reduc- 
tion, is an interesting way to control the iteration process and to help identify convergence 
conditions for the perturbation expansion. 



IV. THE RESULT TO SECOND ORDER 

By applying the procedure introduced in last section, the second order generator of the 
minimal guiding-center transformation is obtained 



(x 2 ) q := 



, n , , „ ab'a — cb'c\ 
b -2c/>cb'b + (22) 



8 

cc — aa 



+ cf> 2 — - — b'b +| (4cab'a + 7acb'a - 9aab'c)] , 
(X 2 ) v := [^(^cb'c-^ab'a + Sacb'a + Saab'c) 

+ 3 {^ T ^b'bb'b} (23) 

, 2 f , _,,20cc + 9aa _, ,7ca - 16ac\ , ., " 
+ 2 < -cb "bb + ( - cb' h ab' ) b'b 

fab"a -cb"c, , ,. ,, x 9ab'a + 7cb'c 
+ 1 | g b + (ab'a - cb'c) 

aa — cc , , , N 3ab'c — 5cb'a 
b'bb'b + cb'a + ab'c )— - — 

8 V ; 32 

+ g{8cb"cc - 5cb"aa + llab"ca 

+ ( - 16cb'cc + lOcb'aa - llab'ca - llab'ac)b'b| 

where is a shorthand for cot tp. To make expressions easier to read, we used the primed 
notation for spatial gradients, and the over-bar over a vector c or a means matrix transpose, 
e.g. ab"ab'c means c • (Vb) ■ [(a ■ V)(Vb)] • a. This notation is close to Littlejohn's notations 
(a be) in 13 , but it is more suited for higher order derivatives. 

In equation (22), the first line is the b-component, parallel to the magnetic field, and the 
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following two lines are perpendicular to b. In equation (23), the first line contains the terms 
depending on V B and the following lines are organized as a polynomial of the pitch-angle, 
or rather its cotangent 0. Note that the terms depending on b ■ VB do not appear in the 
first line, because they are rewritten using formula (21). The expressions involved in such a 
second order result can be written in many different but equivalent ways, as already noticed 
by Northrop and Rome, and a rule for standardizing them is needed for the derivation to 
be efficient. 

If there are many terms, this is especially because each term appears several times, with 
permutations of a and c ; these permutations are often condensed into one single tensor. For 
instance, the last line can be written just 

Vb : n- (b'b), 

by defining the triadic tensor 

II := — 16ccc + lOaca — llcaa — llaac 

Then the number of terms is strongly reduced, each order in <p has one or two terms, which 
shows that the result is not that complicated actually. However, we chose to avoid intro- 
ducing such intermediate quantities, since they make formulae shorter but less explicit. 



The second order vector field gives the minimal guiding-center change of coordinates to 
second order 

[l-(L 1 )-(L 2 )+(^)+0(6 3 )]z, 

where the vector fields Xi and X 2 generating the Lie-transforms l_i and L 2 are given by 
formulae (15), (16), (22) and (23). 

So, to second order, the Larmor radius is 

(xOz-^XiV 



r i =(z-z) q = (X 1 ) q + (X 2 ) q - 



+ 0(e 



p sin t£ . 
eB ' 



(24) 



b -20cb'b + 



5ab'a — cb'c 

8 



+ af^ + 2 ^ + 0a(cb'a-ab'c) 
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And to second order, the reduced pitch-angle is 



<P 



1-(U)-(L 2 



0( t 



_ p sin <p 
~r e B 



-0ab'b - 



cb'a + ab'c 

4 



(25) 



-( £ S £ ) 2 [ilB{- 6 ^b , b + 2ccb'c 

— 2cab'a + acb'a + aab'cj 



+ 



3 r aa cc , 



'{ 



b'bb'b} 



<P 2 



-cb"bb + -5cb 



,2cc + aa 



+ 3ab 



,ca — 2ac 



b'b 



^5aWa-cWc + ga + « Wb 

8 8 

3(ab'a - cb'c) 2 . , . x ab'c-3cb'a 

— + (cb'a + ab'c)- 



32 
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+ g{4cb"cc-cb"aa + 7ab"ca 

+ ( - 8cb'cc + 5cb'aa - 4ab'ca - 7ab'ac)b'b|] . 

The ordering in the Larmor radius is obvious. The first line corresponds to the zeroth and 
first order terms, and all the following lines are the second order terms, which are organized 
as in formula (23). 

We wrote all these formulae exactly as they are yielded by the procedure, because it 
illustrates both the mechanism of the derivation and the structure of the resulting change 
of coordinates. 

Indeed, all the above formulae rely only on a very restricted alphabet of entities : B, 
b, c, a, V and the variable ip, or more precisely (j), if we discard the simp occurring in the 
pre-factor (Larmor radius) and the quantities p and e, which are mute parameters in the 
derivation. 

Each formula is a series in the Larmor radius tl — or more precisely in the ordering 

parameter of guiding-center reduction e = r^V, i.e. the magnetic inhomogeneity at the 
scale of the Larmor radius. This dimensionless parameter is not just a number, it is to 
be understood in the sense that the only dimensional quantities involved are the Larmor 
radius, appearing as a pre-factor to the power given by the expansion order, and the gradients 
(acting on the magnetic field B and b), with the same order as the expansion order. 
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At each order in e, formulae are polynomials in each variables (except B whose disposition 
obeys trivial dimensional rules). Using this fact makes the derivation much simpler. Finally 
the iteration at each order consists in one single formula (8) with two elementary operators 
(derivatives and gyro-integral) ; in addition, they are applied on terms composed of very 
few elementary entities, and those terms are just polynomials. 

As the order in the small parameter e grows up, the polynomial increases its order in each 
of the variables. From the iteration mechanism, a rough estimate shows that the n-th order 
generator X n should be a polynomial of order at most 2n_1 b 3n_1 c 2n , in addition to being a 
monomial in r£V™. Such features in the reduction transformation are useful to consider when 
the asymptotic behaviour is addressed. For instance, a hypothetical convergence condition 
would clearly involve the expected condition on the magnetic inhomogeneity at the scale of 
the Larmor radius, but the role of suggests conditions on the pitch-angle as well. This 
means that in the guiding-center reduction, the direction of the particle velocity must not 
be too close to the direction of the magnetic field. This is in complete agreement with 
the physical intuition, but the polynomial behaviour can help make this intuitive statement 
more precise. 

Another interesting feature is well emphasized by the formulae as they are written : there 
is a link of parity between the order in 0, the order in tl and the fast-angle harmonic (which 
corresponds to the order in c or a). All non-zero terms of order r % L (fP are harmonics of parity 
(— in c for X q , and harmonics of parity (— in c or a for X^. This parity relation 
could already be observed in the first order results (15)-(16), although it was not so obvious, 
because few terms were present. 

It seems that the polynomiality in and the parity relation were first noticed in 23 , where 
they were observed in the derivation of the magnetic moment series, and originated from the 
structure of the operator to be inverted for the secular equation. Here, they are observed 
in the guiding-center reduction, and appear as originating directly from the equations of 
motion (3) or (13), and from the action of derivatives, especially <9^0 and Vc. 

The second order results can be expressed using well-known quantities in guiding-center 
works. For instance, in the term of order 1 in formula equation (25), cb'a + ab'c and 
cb'c — ab'a can be recognized, commonly written — 2Vb : a.i and 4Vb : a 2 where the dyadic 
tensors ai = — ^4^- and a 2 = are well-known in guiding-center works 1 ' 7 , and they were 
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already met in first order results (16). 

Also, in the last line of equation (25), each of the terms is harmonic of order 3 in the fast 
angle c or a. The sum can be rewritten 

( - 3r m a + ab'(bx ) - 8V • b c + 2cb'a a) b'b , (26) 

where three out of four terms have been combined to obtain an harmonic of order 1 in the 
gyro-angle, which is the maximum we can do. Then, in the coefficients, the divergence of 
magnetic field lines V • b = cb'c + ab'a is recovered, as well as their twist 

T rn := b • V x b = (b x V) • b = ab'c - cb'a , 

which quantities are most used in guiding-center theory. 

Such an expression as (26) is useful to get a physical intuition of the terms involved, but 
it does not emphasize as much the mechanism of the derivation nor the polynomial structure 
of the results. Furthermore, it is not unique, e.g. it could equally be written 

( - 5r m a - ab'(bx) - 6V • b c - 2cb'c c)b'b, 

or even 

( - 13r m a - 9ab'(bx) - 8cb' + 2ab'ac) b'b , 

or (at least) nine other equivalent ways of writing this expression. In all of them, three out 
of four terms have been combined to obtain harmonics of order 1 in the gyro-angle, and only 
the last term remains an harmonic of order 3. 

So, we preferred to give the expressions as they come out from the procedure, decomposed 
in elementary terms in a unique, standardized way, as is needed for an algorithmic procedure, 
especially when it can generate very many terms. 

Anyway, it is quite common in perturbation theory and well known in guiding-center the- 
ory that when the order increases, formulae become messy and complicated, and the physical 
interpretation of each term loses some of its relevance 6 . The main point is to standardize 
the derivation and the involved expressions, to make things tractable and as clear as possible. 

Compared to the literature, our fluctuating first order results are exactly identical to 
the usual results of guiding-center derivations, e.g. in 7 . This is because this part of the 
reduction is unique for all guiding-center reductions, since the arbitrariness is only in the 
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average transformation generator. Our second order results are almost identical to previous 
results. For instance, the Larmor radius (24) is exactly identical to the standard result of 7 , 
except for the term <fi 2 , which was absent from Littlejohn's result. As for the pitch-angle, 
the first order term in (25) is exactly the same as in 7 , except the averaged term, which is 
absent from our result. 

A difference with previous results is not a surprise because here the minimal guiding- 
center reduction is considered, and the averaged part of the transformation generator X have 
been set to zero. It is just an effect of the guiding-center reduction non-uniqueness, which 
was well emphasized in 6 and 16 : even between the classical derivations by Lie-transforming 
the phase-space Lagrangian, such as 7,13 ' 15,16 , various changes can be identified in the choices 
adopted for the second order reduction. 

In previous publications, the second order pitch-angle generator (X 2 ) ¥ > was not computed 1 ' 6,7 , 
because guiding-center reductions by Lie-transforming the phase-space Lagrangian requires 
quite a lot of algebra to get (X 2 ) ¥ ,. Even the recent work 15 explicitly chose not to compute 
it, even though it aimed at improving the second order terms and already made quite a lot 
of computations. Finally, (X 2 ) ¥ > was obtained in the soon to be published work 16 , which uses 
a special property between the second and first order terms to avoid pushing the derivation 
to second order. 

These difficulties can be partially explained because when Lie-transforming the La- 
grangian, there is some mixing between the orders in e of the various quantities 1,13 . Thus, 
some components of the first order transformation generator X x are identified at the first 
order analysis of the reduced Lagrangian f i ; other components of Xi are determined at 
second order analysis F 2 , at the same time as some of the components of X 2 ; and the last 
components of Xi are determined at third order T3, at the same time as some components 
of X 2 and of X3. It clearly makes the scheme more involved. This is very different from 
what happens when Lie-transforming the equations of motion, where there is no mixing 
between the order in e of the various quantities, and the non-trivial expression for (X 2 ) v 
was obtained by a direct derivation to order two with very limited algebra. 

Let us turn now to the guiding-center equations of motion. Equations (4) and (5) show 
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that to second order, the drift equations are given by 



e L z 



K 



K (e L z) 
1 + (L 1 ) + (L 2 ) + 



+ 0(e 3 



(27) 



This is just the zero Fourier component of e L z, which was put to zero when computing 
the action of osc in (4). So, this term was already computed in the derivation. 

In this sense as well, Lie-transforming the equations of motion is more straightforward 
because one actually directly derives the reduced motion, whereas when Lie-transforming 
the Lagrangian, one derives the reduced Lagrangian ; the reduced equations of motion must 
be obtained in a second step, for instance by formula (27), which involves some algebra 
because of the lengthy expressions at order 2. 

So, the reduced equations of motion are obtained here with no additional computation 

as 



(28) 



+b»^ + ^(ac-ca)b / b)] 



_ p sin tp 
m 



b(p + r L 



2 

bxVB 
2B 



>Vxb 



+ 2 b x b'b 



[- 



cb'c+cb'a 
2 



(29) 



B' ac— ca l./ 



b'b + 



a(b'b)'c-c(b 



_ p sin ip 



eB r 2 

-^r + r L (j) f^b x b'b + 



b^b)j 



In these equations, we first gave the expression provided by the algorithm and then 
rewrote it using physical quantity, in order to show that the equations are indeed gyro- 
averages, and do not depend on the fast variables c and a. 

The right-hand side of equations (28)- (29) is expressed in the reduced variables (e.g. 
the field is evaluated at the guiding-center position q and the pitch-angle is actually the 
reduced pitch-angle 0), but we dropped the bars for simplicity. Notice that the right-hand 
side of equations (28)-(29) is expressed in the initial variables. It is why equations (22)-(25) 
indeed define a coordinate transformation, whereas equations (28)-(29) indeed define the 
guiding-center dynamics 7 . 

In results (28)-(29), the lowest order term is again exactly identical to the usual guiding- 
center reductions, and the second order drifts are almost identical. The guiding-center 
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dynamics (28) contains exactly the four expected terms : the parallel motion along the 
magnetic field lines, the grad-B drift, the curvature drift, and the Banos drift. But the 
Bahos term b b V 2 xb was not in the reduced dynamics of 7 . This is again because the minimal 
reduction does not aim at minimizing the average dynamics, but at minimizing the change 
of variables. To make it short, it can be said that the reduction of 7 puts the Banos term into 
the change of variables X whereas the minimal reduction lets it in the equations of motion. 
This difference is further considered in the next section. 



V. ON NON-MINIMAL GUIDING-CENTER REDUCTIONS 

In case there is an electric field E in addition to the strong magnetic field, the initial 
equations of motion (13) become 



z = 



eB 



+ 



p sin <p 
eB 



b cot ip+c 


-(b cot y>+c)-Vb-c 




E-psinip 
+ Bp simp | E-(cos V E-b) 
(E-a)a 



As for the variable 6, its dynamics becomes 



e 



eB 



+ cot <p R ■ Vb • a + • Vei • e 2 - 

~ 777 777, 



eE-a 

p simp 



In these time evolutions, the two common small parameters of guiding-center theory 
appear : in strong magnetic field, only one term implies fast dynamics, it is of order ul 
among the other terms, the ones that do not depend on the electric field are of order 
and the ones that do depend on E are of order e si ^ . Thus, two ordering parameters are 
now involved, the magnetic parameter e := = ^ is the same as in the case of a pure 



magnetic field ; the second parameter e# 



uiLm eB 
e E 



is induced by the electric field. 



u>lp simp v±B 

It is the electric force over the magnetic force, and is assumed to be small in guiding-center 
theories 1 ' 2,15 . 

Under this assumption, the procedure for the minimal guiding-center reduction is un- 
changed ; the presence of an electric field only causes some additional contributions in the 
term N n of equation (8), and makes the transformation include one more coordinate : as 
the norm of the particle momentum is no more conserved, this coordinate must be changed 
exactly in the same way as the coordinates (q, <p) in the previous sections. Notice that 
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this minimal guiding-center reduction will imply p to have slow dynamics, but not to be a 
constant of motion. We turn to this point in a few lines. 



In the previous sections, the minimal guiding-center reduction appeared as a natural 
choice when Lie-transforming the equations of motion. Other possibilities are available. 

For instance, one can think of treating the gyro-angle 9 as the other coordinates, and 
changing coordinate 9 — > 9 to gyro-average its dynamics as well. Thus, all of the reduced 
dynamics would involve only slow variables. Especially, the reduced gyro-angle 9 would have 
fast dynamics, but would depend only on the slow dynamics, since it would be independent 
of 9 itself : 9 = #(q, ip,p). So, once the motion of the slow variables (q(£), is known, 
the fast dynamics would be trivial to integrate : 

e(t) = e(to) + J*dt§(q(t),m,p)- 

The corresponding requirement relies on equation (4) written for the gyro-angle compo- 
nent 

= osc (e Ll+L2+ -z) e , (30) 

which is to be solved for L e (or rather X#). The operator to be inverted is exactly the same 
as in the previous sections, since 

-osc(L n+ i(z_ 1 )) 

= -osc^(X n+1 ) i a i (z_ 1 ) e - (z_ 1 ) e a e (X n+1 )^ 

= f d e (X n+1 ) + o.t. , (31) 

where o.t. means other terms that are already known : they do not contain the unknown 
(X. n+ i) e because (i-i)e does not depend on 9. So, just as in the previous sections, a solution 
of equation (30) can be found at arbitrary order in the small parameter, and the dynamics 
of 9 can be averaged as well as the dynamics of (q, tp), as is generally done in guiding-center 
reductions. 

Notice that the averaging process for the reduced coordinates (q, </?) is possible without 
having to deal with resonances 22 because there is no small divisor since the fast part of the 
dynamics involves only one angle 9, as illustrated by formula (8). This is a rather broad 
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property. Now, for the dynamics of 9, the averaging process is possible only because the 
lowest order term in the motion does not depend on 6 itself, as is illustrated by formula 
(31). This is a more specific property. 

For instance, after the guiding-center reduction has averaged the motion over the gyro- 
angle, the bounce- angle can be used to bounce-average the motion 1 . For this reduction, the 
procedure of section 2 can be applied to remove the bounce time-scale from the dynamics 
of the three other components, but averaging the dynamics of the bounce-angle coordinate 
6h as well is guaranteed only if the lowest order bounce-angle dynamics is independent of 
the bounce-angle, which is a condition over the definition of 9 h . 

A second natural requirement concerns the average part of the change of variables avg(X). 
In the minimal derivation, it is trivially put to zero. In the case of a non minimal reduction, 
it appears as an additional freedom in the reduction process, which allows for additional 
requirements, to be chosen. 

For instance, the transformation being a pure gyro-fluctuation avg(z — z) = would be 
interesting : once the effects of the fast coordinate averaged, then the transformation would 
become zero ; so, the reduced coordinates would remain in some sense close to the initial 
ones. 

In the minimal guiding-center reduction, this result is almost obtained, since the trans- 
formation generator X is already a pure gyro-fluctuation, but the whole transformation is 
not a pure fluctuation, because the non-linear terms in X have non-zero gyro-averages (see, 
e.g. formulae (24)- (25)). 

This can be corrected by defining a non- zero average avg(X) to cancel the average of the 
non-linear terms. Indeed, the requirement 




implies 



avgl Lz j = avg (e L z — (1 — L)z) . 
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Now, in the equation at order n 

avg(X„) = avg(Yz) = avg(V L z - (1 - L)zj 

= avgff«411z) (32) 

= avg ^ U- lLl+ U_ 2 L 2+ ... z _ L„_ 2 L 1 L 1 + ... Z + \ 

the right-hand side depends only on generators of orders lower than n. So formula (32) is no 
more an equation, but just a definition of avg(X). Up to second order, equation (32) writes 

avg(X 1 ):=0, 
avg(X 2 ) := avg (k±z"j . 

The right-hand side was already computed in the minimal derivation, so that the additional 
requirement does not imply an additional computation. 



Another possible requirement, commonly considered, is to include the magnetic moment 
/i among the reduced coordinates as a conserved quantity instead of p. Indeed, the basic 
conserved quantity p is not an adiabatic invariant in general when an electric field is present, 
whereas n is so. The requirement is to put to zero not only the fluctuating part of the 
reduced motion osc(z)^, as in equation (4), but its averaged part as well avg(z)^ = 0, which 
simplifies the reduced dynamics in a drastic way. 

For the other coordinates, the same requirement can not be asked, because there is no 
constant of motion independent of p and /x, since otherwise the Hamiltonian motion with 3 
degrees of freedom and 3 independent constants of motion would be integrable. For other 
coordinates, one can not ask so strong a simplification of the reduced dynamics avg(z) for 
all orders in e, but one can ask it for the orders higher than two for instance. In this way, 
the reduced motion would be exactly known to all orders even before computing the change 
of variables. Of course, this is not guaranteed, since it is possible only if the chosen reduced 
dynamics is equivalent to the initial particle dynamics 22 . 

As an example, the minimal guiding-center reduction lets an average contribution (29) 
in the equations of motion for the reduced pitch-angle. Perhaps a stronger reduction could 
cancel it at least for the orders higher than two, by choosing a convenient (non-zero) averaged 
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part X for the transformation generator. More precisely, in the equation at lowest order, 
the unknown X x has no influence on z , which means that (p can not be made a constant 
of motion, as was expected. But at the following order, X_i now contributes, since at that 



where we used that avg (L X2 z_i) v = ; X and X denote respectively the averaged and the 
fluctuating part of X ; and the symbol o.t. is a shorthand for all the other terms, which do 
not contain X x . 

Equation (33) shows that, with the additional requirements, the operator to be inverted 
for X can be quadratic. Even when it is not, e.g. at the next order, the linear operator may 
not be trivial ; a part of it will be given by L^Zq = (X. l )jdjZ — (z )jdj'X 1 , whose coefficients 
have the non-trivial expression given by (3). Such an operator may not be invertible : if 
the right-hand side is not in its range. This is the classical problem of secular terms in 
perturbation theory. Even when the operator is invertible, an explicit inverse may not be 
obvious to get. 

Less strong requirements can be considered than putting to zero the higher order reduced 
dynamics. For instance, Northrop and Rome mention the possible requirement b • q = 0, 
i.e. the guiding-center motion is a pure drift across the magnetic field lines. In any case, 
care must be taken about additional requirements. They imply differential equations that 
are not so simple to deal with, and that even can be impossible to solve. It reminds as 
well that the minimal guiding-center reduction is very nice with its trivial operator or 
equivalently — • d c . 

The requirement aiming at the magnetic moment plays a special role. Its existence can 
be viewed as a consequence of the averaging reduction, as shown by Kruskal 4 . The corre- 
sponding secular differential equation can be solved at any order, but it is outside the scope 
of the present paper and will be the topic of another paper 23 . 

In the approach based on Lie-transforming the equations of motion, the additional re- 
quirements do not change the equation (4) that gives the fluctuating components of the 
transformation ; they generally aim at a further simplification of the reduced dynamics and 
imply additional differential equations such as (33), which impose the averaged components 



order, the equation = avg (e Ll+L2+ -"z) writes 




= o.t. 



(33) 
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of the transformation. Those equations are not easily solved, and from this point of view, 
the method of Lie-transforming the phase-space Lagrangian is more powerful, because it 
does not rely on differential equations, but on algebraic equations, easier to study. 

For instance, the idea to put to zero the averaged components of the reduced motion (at 
least at orders higher than 2 or 3 in the ordering parameter) is included in usual guiding- 
center derivations, which obtain that it is possible for six components of the Lagrangian, out 
of seven. This kind of result would not be so easy to get by Lie-transforming the equations 
of motion. 

Another additional requirement regards the Hamiltonian structure of guiding-center dy- 
namics. The initial motion z is Hamiltonian, and the reduced motion should be Hamiltonian. 
Indeed, the full reduced motion z is Hamiltonian, even if it has a different Poisson bracket 
as the initial motion, since it is just given by a (non-canonical) change of coordinates. 

However the true reduced motion involves only the slow variables(q, <p), and the corre- 
sponding dynamical system (q, 0) is not guaranteed to be Hamiltonian, because it is given by 
a truncation of the full reduced dynamics, and truncations do not preserve the Hamiltonian 
structure in general. 

The preservation of the Hamiltonian structure for the 4-components guiding-center dy- 
namics can be considered as an additional requirement. It is hard to deal by Lie-transforming 
the equations of motion, because when deriving reduced models by working on the equations 
of motion, the Hamiltonian character is not worked on, it is observed a posteriori as preserved 
or not. On the contrary, a Hamiltonian 4-components guiding-center dynamics is easily ob- 
tained when working on the phase-space Lagrangian, for instance by Lie-transforming it 
in such a way that the reduced Poisson bracket is quarter-canonical. This is what most 
guiding-center reductions do. 

VI. CONCLUSION 

Lie-transforming the equations of motion gives a systematic procedure for the guiding- 
center reduction at any order in the Larmor radius expansion. The procedure is very effi- 
cient : it just writes the minimal requirements for the reduction, expands the equation in the 
Larmor radius, and inverts the generator of Larmor gyration to get both the change of coor- 
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dinates and the reduced equations of motion order by order. The full second order reduction 
was straightforwardly obtained, in contrast with previous derivations of the guiding-center 
reduction. 

The corresponding transformation generator is a pure gyro-fluctuation and only four of 
its components are non-zero. This is the bare minimum since exactly the fluctuating part 
of the slow reduced motion involves a fast time-scale that has to be removed. In this sense, 
the reduction is minimal and unique. All the arbitrary components of the transformation 
generator are set to zero. 

The results bring insights into the structure of the guiding-center formulae, which show 
up a polynomiality in b, c, a, r L , V and <fi, and a parity relation between the orders in r L , <f> 
and the fast angle c or a. This polynomiality makes easier the algorithm, which consists in 
applying at each order two operations (derivation, and gyro-integration) onto a polynomial 
of a very restricted alphabet of entities. 

The method can be applied to perform the bounce- average reduction. Also, when an 
additional electric field is present, the procedure is exactly unchanged even if formulae have 
additional terms. 

For the gyro-angle, no gyro-gauge was introduced, since the initial physical coordinate c 
was used. This could bring interesting contributions to some issues involved in the traditional 
gyro-angle 9. Indeed, a recent work 21 was lead to a similar orientation when studying the 
questions raised by the non-existence of a global gauge to define a gyro-angle. 

Thus, this approach of the guiding-center reduction, with its minimal procedure and 
results, can contribute to a better understanding of the guiding-center reduction by giving 
a simplified point of view. 

Non-minimal reductions can be considered by imposing additional requirements for the 
guiding-center reduction. 

The magnetic moment is commonly included in the reduced coordinates. This can be 
done when working on the equations of motion by solving a secular differential equation, 
but it introduces new features and will be reported elsewhere 23 . 

Some averaged terms can be transferred from the reduced dynamics to the change of 
variables, in order to get a stronger guiding-center reduction. This means having non-zero 
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gyro-averaged components in the transformation generator, and these components are again 
obtained by solving secular differential equations. From this point of view, Lie-transforming 
the phase-space Lagrangian is often easier, especially because it mainly relies on algebraic 
equations. Also, it naturally induces a Hamiltonian structure for the 4-components guiding- 
center dynamics. 
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